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(Dated: March 7, 2008) 

The momentum distribution in a Fermi gas with two spin states and a large scattering length has 
a tail that falls off like 1/k 4 at large momentum k, as pointed out by Shina Tan. He used novel 
methods to derive exact relations between the coefficient of the tail in the momentum distribution 
and various other properties of the system. We present simple derivations of these relations using 
the operator product expansion for quantum fields. We identify the coefficient as the integral over 
space of the expectation value of a local operator that measures the density of pairs. 
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Many-body systems of fermions have long been of great 
importance in astrophysics, nuclear physics, and solid 
state physics. The development of trapping and cooling 
techniques for ultracold atoms has made them important 
in atomic physics as well. In this case, the strength of the 
interaction is governed by the 2-body scattering length 
which can be controlled experimentally, adding a new 
dimension to the problem 

If the scattering length a is much larger than the range 
of the interactions, the system has universal properties 
that are determined only by the large scattering length. 
For sufficiently low number density n, the universal prop- 
erties can be calculated using perturbative methods. If 
n\a\ 3 is comparable to 1 or larger, the problem becomes 
nonperturbative. In the special case of two equally- 
populated spin states, systematically improvable calcu- 
lations are possible using Monte Carlo methods. If the 
populations are not equal, this approach suffers from the 
fcrmion sign problem. If there are 3 or more spin states, 
the problem is complicated by the Efimov effect Q . The 
challenging nature of the general problem makes exact re- 
sults very valuable. One case in which exact results are 
known is the unitary limit a = ±oo, where they can be 
derived by exploiting scale invariance Q and conformal 
invariance 

In 2005, Shina Tan pointed out that the momentum 
distribution in an arbitrary system consisting of fermions 
in two spin states with a large scattering length has 
a large-momentum tail that falls off as 1/k 4 [5j. The 
number of fermions with momentum larger than K ap- 
proaches C/(ir 2 K) as K — ► oo , where C depends on the 
state of the system. Tan used novel methods involving 
generalized functions to derive exact relations between C 
and several other properties of the system. An example 
is the adiabatic relation that gives the change in the total 
energy E due to a small change in a 0: 



dE 

d(l/a) 



h 2 

Airm 



C 



(1) 



Tan referred to C as the integrated contact intensity, 
which we will abbreviate to contact. The Tan relations 
hold for any state of the system: few-body or many- 
body, homogeneous or in a trapping potential, superfluid 
or normal, zero or nonzero temperature. 

In this Letter, we show that the Tan relations can be 
derived using the operator product expansion (OPE) for 
quantum fields. The OPE was proposed by Ken Wilson 
in 1969 as a formalism for dealing with the strong in- 
teractions associated with the nuclear force. The OPE 
has become a standard tool to understand the relativistic 
quantum field theories that describe elementary particles 
We will apply the OPE to the strongly-interacting 
nonrelativistic system consisting of fermions with two 
spin states and a large scattering length. We identify 
the contact C as the expectation value of the integral 
over space of a local operator that measures the density 
of pairs. The OPE provides new insights into the Tan re- 
lations and makes it easier to generalize them to systems 
with additional degrees of freedom or more complicated 
interactions. 

A system consisting of fermions in two spin states la- 
belled by a = 1,2 can be described by a quantum field 
theory with two quantum fields vb a (r). The number op- 
erator is J2a I d 3 R ^J.^ct(-R) \M- The momentum distri- 
bution p a (k) for fermions with spin a can be expressed 
as 

p a (k) = [d 3 R[d 3 re ihr (iPt(R-±r)i> a (R+hr)) . (2) 

Its behavior at large k is determined by the matrix el- 
ement at small r. If the fermions are non-interacting, 
the quantum fields can be expanded as Taylor series in 
r. The resulting expansion for p a (k) can be expressed in 
terms of the Dirac delta function in k and derivatives of 
the delta function. This indicates that p a {k) has a finite 
range in k. Simple examples are an ideal gas of fermions 
at temperature, for which p a {k) vanishes if |fe| is larger 
than the Fermi momentum, and an ideal gas of fermions 






FIG. 1: (Color online) Integral equation for the scattering of 
a pair of fermions. The blob represents the amplitude iA(E). 
The solid dot represents the interaction vertex —ig(A). 



at high temperature, for which p a {k) is a Gaussian func- 
tion of k. 

If there are interactions between the fermions, the ma- 
trix element in Eq. ([2]) may not be an analytic function 
of r at r = 0. The momentum distribution p a (k) may 
therefore have a large-momentum tail that falls off like 
a power of k = \k\. Tan showed that for fermions with 
two spin states and a large scattering length, the tail is 
proportional to 1/fc 4 and is the same for both spin states: 



p a {k) — > C/fc 4 



(3) 



where C is the contact. This power-law behavior can 
arise from a term in the matrix element in Eq. ^ that 
is linear in r = \r\ as r — > 0. It holds for all states that 
include at least one particle of each spin. The dependence 
on the state enters only through C. 

The possibility of a power-law tail in the momentum 
distribution can be understood from Wilson's OPE. Wil- 
son proposed that in a quantum field theory a product 
of local operators separated by a short distance can be 
expanded in terms of local operators. The expansion of 
the product of quantum fields in Eq. @ has the form 



(4) 



where the sum is over local operators O n (R) that can be 
constructed out of the quantum fields and their gradients. 
They include the operators that arise from multiplying 
the Taylor expansions of the operators on the left side of 
the equation, but there can be other operators as well. 
The functions C^nfV) are called Wilson coefficients or 
short- distance coefficients. Some of them can be nonan- 
alytic at r = 0, so the momentum distribution p a {k) can 
have a power-law tail. 

The Hamiltonian density for the quantum field theory 
that describes fermions with two spin states and a large 
scattering length in an external potential V(R) is 



(5) 



where V = V(R) J2a VaVv- For simplicity, we have set 
h = 1. The superscripts (A) on the operators in Eq. ([5|) 
indicate that their matrix elements are ultraviolet diver- 
gent and an ultraviolet cutoff is required to make them 
well defined. For the ultraviolet cutoff, we impose an 



(a) 




(b) 




(c) 




((/) 




FIG. 2: (Color online) Diagrams for matrix elements 
of the operators (a) ipl(— \r)ij) a {+\r), (b) ^J^v(O), (c) 

^/)J'0|VlV'2 (0)> an d (d) ip\ipi{—^r)iplip2(+^r) between scat- 
tering states consisting of a pair of fermions. The open dots 
represent the operators. For each diagram, there are three 
other diagrams in which there is no scattering of the two in- 
coming fermions or the two outgoing fermions or both. 



upper limit |fc| < A on momentum integrals. In the 
limit A — > oo, the Hamiltonian density in Eq. ([5]) de- 
scribes fermions with zero-range interactions and scat- 
tering length a if we take the coupling constant to be 



5(A) = 
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The amplitude for the scattering of a pair of fermions can 
be calculated by solving the Lippmann-Schwinger inte- 
gral equation, which is represented diagrammatically in 
Fig. [1] The solution A(E) depends on the total energy 
E of the pair of fermions in the center-of-mass frame 
and not separately on their momenta. After substituting 
Eq. ([6]) for g (A), the solution in the limit A — > oo is 



A{E) = 



■1/a + ^J-mE ■ 



(7) 



The T-matrix element for scattering of a pair of fermions 
with momenta +p and —p is obtained by setting E = 
p 2 /to. 

We now proceed to show that the OPE in Eq. Q in- 
cludes a term with the operator O12 — ipliplipi^^ an d 
that its Wilson coefficient 6^12 (r) is linear in r = |r|. 
Since the OPE is an operator equation, C a> v2.{r) can 
be determined by calculating the expectation value of 
both sides of the OPE in the simplest states for which 
(ipliplipiip^) is nonzero. We choose the state | ± p) 
consisting of two fermions with spins 1 and 2 and with 
momenta +p and —p. 

We first consider the left side of the OPE in Eq. ((4}. 
The operator product at R = can be represented dia- 
grammatically by a pair of dots that correspond to the 
points where it annihilates an incoming fermion and 
— \r where it creates an outgoing fermion. The expecta- 
tion value can be represented by the sum of the diagram 
in Fig. Ufa) and the three diagrams with one or no scat- 
terings. The contributions from the other three diagrams 
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are analytic at r = 0. The contribution from the diagram 
in Fig. [2a) is 

{^{-\r)M^r)) ±p ^ ) =im 2 A 2 { P 2 /m)e^/{^p) . 

(8) 

If we expand this in powers of r, the terms with odd 
powers of r are nonanalytic at r = 0. 

We now consider the right side of the OPE in Eq. Q . 
The expectation value of ip^ipa can be represented by the 
sum of the diagram in Fig. [^b) and the three diagrams 
with one or no scatterings. The contribution from the 
diagram in Fig. [2jb) is 

(iplM0))±P^ h) = im 2 A 2 (p 2 /m)/(8np) . (9) 

This matches the r° term in the expansion of Eq. © 
in powers of r. Thus we find that the Wilson coeffi- 
cient of ipl.'4! a (R) in Eq. is simply 1 in accord with 
naive expectations. To match the term linear in r in the 
expansion of Eq. 1(5)), we must find an operator whose 
expectation value in the state | ± p) has the momentum 
dependence A 2 (p 2 /m). One can deduce this operator by 
realizing that the nonanalytic behavior at r = must 
arise from the region of the diagram in Fig. Ula) i n which 
there is large momentum flowing in the lines that con- 
nect the operators to the amplitudes and the line that 
connects the two amplitudes. If we shrink all these lines 
to a point, they reduce to a vertex with two incoming 
lines and two outgoing lines. The simplest corresponding 
operator is O12 = V , iV4V , iV'2 • We can verify that this 
is the correct operator by calculating its matrix element, 
which can be represented by the sum of the diagram in 
Fig. (He) and the three diagrams with one or no scatter- 
ings. By using the integral equation in Fig. [I] the sum of 
the four diagrams can be expressed in the simple form 

(V4v4V>i4 A) (0))±p = m 2 g- 2 (A)A 2 (p 2 /m) . (10) 

This has the same dependence on p as the term linear 
in r in the expansion of Eq. {§]), To match this term, 
the Wilson coefficient must be C a 12(f) = — .9 2 (A) r j (8tt). 
The expectation value of this term in OPE can be written 

C ail2 (r)(0 12 (R)) = -L-{g 2 ^l^ 2 (R)) . (11) 

We have attached the factor <? 2 (A) to the operator 

V'lV'IV'iV'a an d suppressed the dependence on A, be- 
cause the resulting operator has finite matrix elements 
in the limit A — > 00, as exemplified by Eq. (|10|) . We 
proceed to use this result to derive the Tan relations. 

Tail of the momentum distribution. The leading 
behavior of Pcr(k) at large fc can be obtained by insert- 
ing the term in Eq. (jlip in place of the matrix element 
in Eq. (|2|). This term gives the asymptotic behavior in 
Eq. with 

C= [d 3 R(g 2 44iPiMR)) ■ (12) 



We will refer to <7 2 V>iV4V'i'02 as the contact density oper- 
ator. This is a positive operator, so C > 0. 

Energy relation. In Ref. [H, Shina Tan derived an 
expression for the energy £ as a linear functional of the 
momentum distributions: 

f d 3 k k 2 f , C\ C 

. . (13) 

This is a functional of p\ and pi 1 because C is de- 
termined by their large-momentum behavior: C — 
limfe— >oo k i p a (k). Tan's energy relation can be derived 
simply by using the expression for g(A) in Eq. ^ to ex- 
press the Hamiltonian density in Eq. © as the sum of 
three terms whose matrix elements are ultraviolet finite: 

a 

+£ 5 2 ^W 2 + V. (14) 

The operator V^J. • V^>^ in the first term has matrix 
elements that diverge linearly as A — > 00. This linear 
divergence is cancelled by the second term, which has an 
explicit factor of A. Integrating over the positions of the 
local operators in Eq. fTi]) , taking the expectation value, 
and using the expression for C in Eq. (|12p . we obtain 
Tan's energy relation in Eq. (fT3)) . 

Adiabatic relation. Eq. (pj can be derived by using 
the Feynman-Hellman theorem: 

dE/da = Jd 3 R(dH/da) . (15) 

Since 7i depends on a only through the coupling constant 
g(A) in Eq. ([6|), its derivative can be written 

dH/da = g 2 i)\^\i> x ^ 2 l{\-Kma 2 ) . (16) 

Upon inserting this into Eq. (fT5|) and using Eq. (fT2)l for 
C, we get the adiabatic relation in Eq. ([T]). 

Virial theorem. The virial theorem for fermions with 
two spin states in a harmonic trapping potential in the 
unitary limit a = ±00 was derived in Ref. [loj ]. Tan has 
derived a generalization of the virial theorem for the case 
of finite a [Til ] . The virial theorem can be derived more 
simply by using the fact that the scattering length a and 
the angular frequency uj of the trapping potential provide 
the only scales for the energy E of a state. Dimensional 
analysis then requires the differential operator (ujd/dio) — 
^(ad/da) to give 1 when acting on Jd 3 R (TL). Using the 
Feynman-Hellman theorem together with Eq. (j 16[) . we 
obtain the virial theorem 

E = 2 Jd 3 R (V) - C/(8nma) . (17) 

We can obtain a simple interpretation of the contact 
density operator by considering the OPE of the number 
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density operators %j}\^>\{— |r) and ^IV^C+s r )- The ex- 
pectation value of their product in the state | ±p) can be 
represented by the sum of the diagram in Fig. Eld) and 
the three diagrams with one or no scatterings. The con- 
tribution that is most singular as r — > comes from the 
diagram in Fig. d^d): m 2 A 2 (p 2 /m)e 2zpr /(I6ir 2 r 2 ). The 
term proportional to r~ 2 has the same dependence on p 
as the matrix element in Eq. (|10|) • Thus the most singular 
term in the OPE is 0] 



V4Vi(#-HV4V>2(i?+|r) 



1 



g 2 44^MR) ■ 



167r 2 r 2 



(18) 



We can define an operator N pail (R,s) that counts the 
number of pairs of fermions with spins 1 and 2 near 
the point R with separation less than s by integrat- 
ing the left side of Eq. (fTB")) over the ball \r\ < s. In 
the absence of interactions, (N pa i r (R, s)) scales as s 3 as 
s — > 0. Eq. (jT8j) implies that in the case of a large scat- 
tering length (iVpair(-R, s)) scales as s 1 . We can interpret 
the contact density operator g 2 ip\iplipiip2 as the limit as 
s^Qof (4ir/s)N pail (R,s). 

One can give a more intuitive interpretation of the con- 
tact C if the fermions with spins 1 and 2 have an inelastic 
two-body scattering channel into other spin states that 
are much lower in energy. In this case, the optical the- 
orem implies that the scattering length a has a nega- 
tive imaginary part. The leading effects of a weakly- 
coupled inelastic channel on low-energy fermions in the 
spin states of interest can be taken into account through 
the small imaginary part of a. The effect on a state with 
definite energy E is to change its time-dependence from 
exp(-iEt/h) to exp(-i(J5 - iT/2)t/%). The probability 
in that state decreases with time at the rate Y/K. The 
adiabatic relation in Eq. (TT|) can be used to derive an ex- 
pression for r to leading order in the imaginary part of 
a: 
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Thus C determines the rate at which low-energy fermions 
are depleted by inelastic collisions. 

In Ref. 0], Tan derived expressions for the contact 
density in the BCS limit (a — * _ ), the unitary limit 
(a — » ±oo), and the BEC limit (a — > + ) by using the 
adiabatic relation in Eq. |T]) as an operational definition 
of C. Our identification of the local contact density op- 
erator g 2 ip\iplipiip2 makes it straightforward to calculate 
the contact density for homogeneous systems directly us- 
ing diagrammatic methods. 

The Tan relations for ultracold fcrmionic atoms with 
two spin states apply equally well to cold neutron matter 



at sufficiently low densities. One advantage of our deriva- 
tion using the OPE is that it makes it easier to generalize 
the Tan relations to more complicated systems. The gen- 
eralization to more complicated interactions is straight- 
forward for any system that can be described by a renor- 
malizable local quantum field theory (l2| . An important 
example is the resonance model that provides a natural 
description of atoms near a Feshbach resonance [l3j . The 
generalization to systems with additional spin degrees of 
freedom, such as nuclear matter at sufficiently low den- 
sities, is complicated by the Efimov effect |2j. In such 
3-body analog of the contact density operator 
may be expected to play an important role. 

The original definition of the contact C in terms of 
the tail of the momentum distribution suggests that the 
Tan relations are relevant only to esoteric aspects of the 
strongly-interacting Fermi gas. However the adiabatic re- 
lation in Eq. {1} and the virial theorem in Eq. (fl7|) make 
it clear that they are actually of central importance. In 
Refs. 0, @, Tan offered suggestions for how these rela- 
tions could be tested experimentally. We leave this as 
a challenge to the ingenuity of experimentalists in cold 
atom physics. 
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